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Measurement of direct photon production at Tevatron fixed target energies 
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Measurements of the production of high transverse momentum direct photons by a 515 GeV/c 7r~ 
beam and 530 and 800 GeV/c proton beams in interactions with beryllium and hydrogen targets are 
presented. The data span the kinematic ranges of 3.5 < pt < 12 GeV/c in transverse momentum and 
1.5 units in rapidity. The inclusive direct-photon cross sections are compared with next-to-leading- 
order perturbative QCD calculations and expectations based on a phenomenological parton-fcy 
model. 

PACS numbers: 13.85.Qk, 12.38.Qk 



I. INTRODUCTION 

Inclusive single-particle production at large transverse 
momentum {pt) has been useful in the development of 
perturbative quantum chromodynamics (PQCD) [1-4]. 
Quantitative comparisons with data have yielded infor- 
mation on the validity of the PQCD description, and on 
parton distribution functions of hadrons (PDF) and the 
fragmentation functions of partons. This paper reports 
precision measurements of the production of direct pho- 
tons with large pt- 

At leading order, only two processes contribute to the 
direct-photon cross section, namely, qq annihilation and 
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FIG. 1: Leading-order diagrams for direct-photon production. 



quark-gluon Compton scattering (Fig. 1). The pho- 
ton's momentum reflects the collision kinematics since 
such photons are produced at the elementary interaction 
vertex. This contrasts with jet production, where the 
hadronization process obscures the measurement of en- 
ergy and direction of the outgoing parton. A complete 
theoretical description of the direct-photon process is of 
special importance as it has long been expected to facili- 
tate the extraction of the gluon distribution of the proton. 
The quark-gluon Compton scattering process shown in 
Fig. 1 provides a major contribution to inclusive direct- 
photon production. The gluon distribution is relatively 
well constrained by deep-inelastic and Drell-Yan data for 
momentum fractions x < 0.1, but less so at larger x [5]. 
Fixed-target direct-photon data were incorporated into 
several previous global parton-distribution analyses [6-8] 
to constrain the gluon distribution at large x, but more 
recent global PDF analyses have not included such pho- 
ton data; this point is revisited in the concluding section 
of this work. 



II. APPARATUS 



A. Meson West spectrometer 
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Fermilab E706 is a fixed-target experiment designed 
to measure the production of direct photons, neutral 
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mesons, and associated particles at high-pr [9-13]. The 
apparatus included a charged particle spectrometer, a 
large liquid argon calorimeter and downstream muon 
identifiers. 

The data collection phase of this experiment spanned 
three fixed-target running periods, including the rel- 
atively low-statistics commissioning run during 1987- 
1988 [14, 15], and the two primary data runs in 1990 
and 1991-1992. This paper reports on data from the 
two primary data runs of the experiment [16]. During 
the 1990 run, the target consisted of two 0.8 mm thick 
copper foils followed by two pieces of beryllium. The up- 
stream piece was 3.7 cm long, while the length of the 
downstream piece was 1.1 cm. In the 1991-1992 run, the 
target consisted of two 0.8 mm thick copper foils immedi- 
ately upstream of a liquid hydrogen target [17], followed 
by a 2.54 cm long beryllium cylinder. The liquid hydro- 
gen was contained in a 15.3 cm long mylar flask, which 
was supported in an evacuated volume with beryllium 
windows at each end (2.5 mm thickness upstream and 
2.8 mm thickness downstream). 

The charged particle spectrometer consisted of silicon 
microstrip detectors (SSDs) in the target region and mul- 
tiwire proportional chambers (PWCs) and straw tube 
drift chambers (STDCs) downstream of a large-aperture 
analysis magnet [12]. Six 3x3 cm^ SSD planes were lo- 
cated upstream of the target region and used to recon- 
struct beam tracks. Two hybrid 5x5 cm^ SSD planes 
(25 fim pitch strips in the central 1 cm, 50 /xm beyond) 
were located downstream of the target region. These 
were followed by eight 5x5 cm'^ SSD planes of 50 fim 
pitch. The analysis dipole imparted a 0.45 GeV/c pt im- 
pulse in the horizontal plane to charged particles. Down- 
stream track segments were measured by means of four 
stations of four views (XYUV) of 2.54 mm pitch PWCs 
and two stations of eight {4XAY) layers of STDCs with 
tube diameters 1.03 cm (upstream station) and 1.59 cm 
(downstream station) [18]. 

Photons were detected in a large, lead and liquid-argon 
sampling electromagnetic calorimeter (EMLAC), located 
9 m downstream of the target [19]. The EMLAC had a 
cylindrical geometry with an inner radius of 20 cm and 
an outer radius of 160 cm. The calorimeter had 33 lon- 
gitudinal cells read out in two sections: an 11 cell front 
section (8.5 radiation lengths) and a 22 cell back section 
(18 radiation lengths). Each longitudinal cell consisted 
of a 2 mm thick lead cathode (the first cathode was con- 
structed of aluminum), a double-sided copper-clad G-10 
radial (R) anode board, a second 2 mm thick lead cath- 
ode, and a double-sided copper-clad G-10 azimuthal (<&) 
anode board. The 2.5 mm gaps between these layers were 
filled with liquid argon. The physical layout is illustrated 
in Fig. 2. 

The EMLAC readout was subdivided azimuthally into 
octants, each consisting of interleaved, finely segmented, 
radial and azimuthal views. This segmentation was real- 
ized by cutting the copper-cladding on the anode boards 
to form either radial or azimuthal strips. Signals from 
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FIG. 2: A drawing of the liquid argon electromagnetic 
calorimeter with some components pulled away in one quad- 
rant to reveal a view of the internal details. 

corresponding strips from all R (or $) anode boards 
in the front (or back) section of a given octant were 
jumpered together. The copper-cladding on the radial 
anode boards was cut into concentric strips centered on 
the nominal beam axis. The width of the strips on the 
first R board was 5.5 mm. The width of the strips on 
the following R boards increased slightly so that the ra- 
dial geometry was projective relative to the target region. 
The azimuthal strips were split at a radius of 40 cm into 
inner and outer segments; each inner strip subtended an 
azimuthal angle of 7r/192 radians, while outer strips cov- 
ered 7r/384 radians. 

The spectrometer also included two other calorimeters: 
a hadronic calorimeter located downstream of the EM- 
LAC within the same cryostat, and a steel and scintillator 
calorimeter positioned further downstream to increase 
coverage in the very forward region. The E672 muon 
spectrometer, consisting of a toroidal magnet, shield- 
ing, scintillator, and proportional wire chambers, was de- 
ployed immediately downstream of the calorimeters [20- 
23]. 

The spectrometer was located at the end of the Me- 
son West beamline. The design of the beamline, primary 
target and primary beam dump were intended to min- 
imize the rate of beam-halo muons incident upon the 
spectrometer. The beamline was capable of transport- 
ing either a primary (800 GeV/c) proton beam or un- 
separated secondary particle beams of either polarity to 
the experimental hall. The beamline Cerenkov detector 
was used to identify the secondary beam particles [24]. 
This 43.4 m long helium-filled counter was located 100 m 
upstream of the experimental target. The positive sec- 
ondary beam with mean momentum of 530 GeV/c was 
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97% protons. The negative secondary beam with mean 
momentum of 515 GeV/c was 99% pions. 

At the end of the beamline was a 4.7 m long stack 
of steel siirrounding the beam pipe and shadowing the 
EMLAC to absorb off- axis hadrons. A water tank was 
placed at the downstream end of this hadron shield to ab- 
sorb low-energy neutrons. Surrounding the hadron shield 
and neutron absorber were walls of scintillation counters 
(VW) to identify penetrating muons. There was one wall 
at the upstream end and two walls at the downstream 
end of the hadron absorber during the 1990 run. An 
additional wall was added to the upstream end of the 
hadron absorber prior to the 1991-1992 run. 



B. Trigger 

The E706 trigger selected interactions yielding high-py 
showers in the EMLAC. The selection process involved 
several stages: beam and interaction definitions, a pre- 
trigger, and high-py trigger requirements [10, 12, 25]. 
A scintillator hodoscope, located 2 m upstream of the 
target region, was used to detect beam particles, and re- 
ject interactions with more than one incident beam parti- 
cle. Additional scintillator with a 1 cm diameter central 
hole was located just downstream of the beam hodoscope, 
and served to reject interactions initiated by particles in 
the beam halo [26]. Two pairs of scintillator counters, 
mounted on the dipole analysis magnet, were used to 
identify interactions in the target. A filter was employed 
for the high-pT triggers to reject interactions that oc- 
curred within 60 ns of one another to minimize potential 
confusion in the EMLAC due to out-of-time interactions. 

For those interactions that satisfied the beam and in- 
teraction requirements, the pT deposited in various re- 
gions of the EMLAC was evaluated by weighting the en- 
ergy signals from the EMLAC i?-channel amplifier fast 
outputs by a factor proportional to sin^j, where 9i was 
the polar angle between the i*^ strip and the nominal 
beam axis. The pretrigger hi requirement for a given 
octant was satisfied when the px detected in either the in- 
ner 128 i? channels or the outer R channels of that octant 
was greater than a threshold value. A pretrigger signal 
was issued only when there was no evidence in that oc- 
tant of substantial noise, significant pt attributable to an 
earlier interaction, or incident beam-halo muon detected 
by the VW. 

Localized trigger groups were formed for each octant 
by clustering the i?-channel fast-outputs into 32 groups 
of 8 channels. Each adjacent pair of 8 channel groups 
formed a group-of-16 strips. If the Pt detected in any 
of these groups-of-16 was above a specified high (or low) 
threshold, then a local hi (or local lo) signal was 
generated for that octant. A single local hi (or sin- 
gle LOCAL LO) trigger was generated if a LOCAL hi (or 
LOCAL Lo) signal was generated in coincidence with the 
PRETRIGGER HI in the same octant. 

Trigger decisions were also made based upon global 



energy depositions within an octant. A global lo sig- 
nal was generated if the total px in an octant exceeded 
a threshold value. The LOCAL0GLOBAL LO trigger re- 
quired a coincidence of the pretrigger hi signal with 
global lo and local lo signals from the same octant. 
The LOCAL LO requirement was included to suppress spu- 
rious global triggers due to coherent noise in the EMLAC. 

The SINGLE LOCAL LO and LOCAL(8)GLOBAL LO trig- 
gers were prescaled to keep them from dominating the 
data sample. Prescaled samples of beam, interaction, 
and pretrigger events were also recorded. 



III. ANALYSIS METHODS 

Data samples contributing to this analysis represent 
an integrated luminosity of 6.8 (1.1) pb~^ for 530 GeV/c 
pBe [pp) interactions, and 6.5 (1.1) pb~^ at 800 GeV/c. 
These samples were accumulated during the 1991-1992 
nm. Results reported in this paper also use 6.1 pb""'^ 
of TT^Bc data recorded during the 1990 run, as well as 
1.4 pb~^ of 7r~Be data and 0.23 pb~^ of -K~p data accu- 
mulated during the 1991-1992 run. The following subsec- 
tions describe the analysis procedures and methods used 
to correct the data for losses due to inefficiencies and se- 
lection biases. Additional details can be found in several 
of our previous papers [10-12, 19]. 



A. Event reconstruction 

Two major aspects of the analysis procedure involved 
charged-particle and calorimeter shower reconstruction. 
The charged-track reconstruction algorithm produced 
track segments upstream of the magnet using informa- 
tion from the SSDs, and downstream of the magnet us- 
ing information from the PWCs and STDCs. These track 
segments were projected to the center of the magnet and 
linked to form the final tracks whose calculated charges 
and momenta were used for the physics analysis. The 
charged track reconstruction is described in more detail 
elsewhere [12]. 

The readout in each EMLAC quadrant consisted of 
four regions: left and right R, and inner and outer 
Strip energies from clusters in each region were fit to the 
shape of an electromagnetic shower determined from de- 
tailed Monte Carlo simulations and isolated-shower data. 
These fits were used to evaluate the positions and ener- 
gies of the peaks in each region. Shower positions and 
energies were obtained by correlating peaks of approxi- 
mately the same energy in the R and $ regions within the 
same half octant. More complex algorithms were used to 
handle configurations with showers spanning multiple re- 
gions. The EMLAC readout was also subdivided longitu- 
dinally into front and back sections. This segmentation 
provided discrimination between showers generated by 
electromagnetically or hadronically interacting particles. 
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FIG. 3: The longitudinal distribution of reconstructed ver- 
tices for 1990 (top) and 1991-1992 (bottom) target configura- 
tions. The data arc represented by the histogram, the results 
from the detailed detector simulation by the points. 



An expanded discussion of the EMLAC reconstruction 
procedures and performance can be found elsewhere [19]. 



B. Selection criteria and corrections 

Only events with reconstructed vertices within the 
fiducial volume of the Be or H2 targets contribute to the 
results presented in this paper. Vertex reconstruction ef- 
ficiencies were evaluated for each target using a detailed 
Monte Carlo simulation of the spectrometer, as described 
in Section III.F below [27]. These efficiencies were used to 
correct for reconstruction losses and resolution smearing 
across fiducial boundaries of the targets. The corrections 
were < 1% for interactions in the bcryUium targets in the 
1990 configuration as well as for the hydrogen and down- 
stream beryllium targets in the 1991-1992 configuration. 
The correction was <3.5% for the upstream Be target in 
the 1991-1992 configuration. The longitudinal distribu- 
tion of reconstructed vertices are shown, for both 1990 
and 1991-1992 target configurations, in Fig. 3. Overlaid 
on the data for the copper, beryllium, and hydrogen tar- 
gets arc the results from the detailed Monte Carlo simu- 
lation. 

Only photons that were reconstructed with the fiducial 

regions of the EMLAC contributed to the cross-section 
results presented. In particular, regions of the detec- 
tor near quadrant boundaries that abutted steel support 
plates, the central beam hole, the outer radius of the 
EMLAC, and octant boundaries were excluded. A sim- 



FIG. 4: Fraction of energy deposited in the front section of 
the EMLAC for identified electromagnetic showers (electrons 

from photon conversions) and hadronic showers (charged pi- 
ous from Ks decays). The data are represented by the his- 
tograms, the results from the detailed detector simulation by 
the points. The distributions have been normalized to unit 
area. 



pie ray-tracing Monte Carlo program was employed to 
determine the correction for the ^-dependent fiducial re- 
quirements. 

Showers that deposited <20% of their energy in the 
front section of the EMLAC were not considered as pho- 
ton candidates since such showers were likely to be due to 
hadronic interactions [19] (Fig. 4). Losses of photons due 
to this requirement were « 2%. A detailed event simula- 
tion (described in Section IILF below) was employed to 
correct for this and other eflfects including reconstruction 
smearing and losses. 

Reconstructed showers were excluded from the photon 
sample when chargcd-particle tracks pointed to within 
1 cm of shower center (Fig. 5). The correction for this 
criterion in the direct-photon analysis is w 1.01 based 
upon studies of the impact of this requirement on recon- 
structed 7r°'s. 

The correction for losses due to the conversion of pho- 
tons into e+e~ pairs was evaluated by projecting each 
reconstructed photon from the event vertex to the recon- 
structed position in the EMLAC. The radiation length of 
material traversed, up to the analysis magnet, was eval- 
uated based upon detailed detector descriptions. The 
photon conversion probability was evaluated and used to 
account for conversion losses. The average correction for 
conversion losses was 1.09 per photon for the Be target in 
the 1990 run (1.08 in 1991-1992 run) and 1.11 per photon 
for the H2 target. 
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FIG. 5: Distribution of the distance squared between the 
reconstructed positions of electromagnetic showers and the 
nearest charged-particle track in data (histogram) and in the 
detailed detector simulation (•). To enhance the electron frac- 
tion in these samples, only showers with > 0.75 of their en- 
ergy deposited in the front section of the EMLAC contribute 
to this plot. The distributions have been normalized to the 
same area. 



C. Trigger corrections 

Trigger corrections were evaluated on an event-by- 
event basis using the energy depositions detected in the 
EMLAC combined with the measured responses of the 
trigger groups in which those energy depositions occurred 
for this particular trigger and event. The probability (pi) 
that the i*'^ trigger group (e.g., a group-of-16 or inner 
128 R channels) in the octant was satisfied was mea- 
sured as a function of the pr reconstructed within the 
group using data samples unbiased with respect to that 
group [25]. The trigger probability for a given octant 
was defined as P = 1 — YI{1 — pi). Trigger corrections 
were evaluated on an cvent-by-cvent basis by applying 
the inverse of the trigger probability as a weight to each 
photon candidate. Photon candidates with P < 0.1 were 
excluded from further consideration to avoid excessively 
large weights. The correction for losses from this require- 
ment was determined from the detailed detector simula- 
tion and absorbed into the reconstruction efficiency de- 
scribed below. 

The results presented in this paper are extracted from 

data selected using several triggers with different thresh- 
olds and prcscalc factors to probe a broad range of pT- 
The lower threshold and prescaled triggers are used to 
populate the lower px range where the data from the 
higher threshold triggers were judged to be less reliable 
(due to larger average trigger corrections) than corre- 
sponding data from the lower threshold triggers. These 
judgements were made by comparing the fully corrected 
results from each trigger in px and rapidity, t/cm- These 
transitions were functions of both pT and Ucm, and 
were different for the three samples. The single lo- 
cal HI trigger was used exclusively for photons with 
PT > 4.0 GeV/c. 



D. Rejection of beam-halo muons 

Spurious triggers were produced by muons in the beam 
halo that deposited significant energy in the EMLAC in 
random coincidence with interactions in the target. Such 
muon- induced showers were interpreted by the trigger as 
high-pT depositions, particularly when they occurred in 
the outer regions of the EMLAC. While the pretrigger 
logic used VW signals to reject events associated with 
such muons, additional requirements on the latched VW 
signals, the total pt imbalance in the event, and the di- 
rection and shape of reconstructed showers were imposed 
during the offline analysis to further suppress this back- 
ground. We note that these criteria successfully elimi- 
nated the beam-halo muon background in the high-pr 
meson samples (see Fig. 6 of Ref. [10] and Fig. 5 of 
Ref. [11]). As expected, the beam-halo muons rates were 
much larger in the secondary beam data samples than in 
the primary beam (800 GeV/c) data sample. 

The Pt balance in the event was defined by calculat- 
ing the Pt of the photons and charged particles which, 
based upon their initial trajectories, would intercept the 
EMLAC in the transverse plane within the 120° sector 
opposite the direct-photon candidate. The ratio of this 
"away-side" pt (P-^™"^) to that of the direct-photon can- 
didate should be near unity in interactions that generate 
a high-pT photon. However, for events triggered by show- 
ers from beam-halo muons, P^^"^ /pt should be near 
zero, since interactions arc typically soft. Showers with 
pojDoy^^^ < 0.3 were excluded from the photon sample. 

A measure of the incident direction of showers ob- 
served in the EMLAC was characterized as 5 = 
-Rfront - (■^front/.Z'back) ^?back, where Pfront and i?back are 
the radial positions of the reconstructed shower in the 
front and back sections, and .^front and .^back are the 
longitudinal locations of those EMLAC sections relative 
to the nominal target location. Showers with 5 ~ were 
generated by particles from the target region, while show- 
ers with large 5 were generated by particles incident from 
other directions. The requirement on 5 was a function 
of radius; showers with S > 0.193 cm for i? < 40 cm 
{6 > 0.0048ii for i? > 40 cm) were excluded from the 
photon sample. 

The electromagnetic shower shape employed in the 
photon reconstruction program was determined for pho- 
tons originating in the target region [19]. Consequently, 
showers resulting from beam-halo muons were frequently 
poorly fit by this shower shape, particularly in the radial 
view. The of the radial view fit therefore provided ad- 
ditional discrimination against beam-halo muons. Show- 
ers with XbJE > 0.1 CeV~^ (where E represents the 
shower energy) were excluded from the photon sample. 

Distributions of P^^'^^ /px, S, and X%JE are shown in 
Fig. 6 for a beam-halo muon enriched sample (left) and 
a photon enriched sample (right). The rejection crite- 
ria discussed above suppressed beam-halo muon-induced 
background while having little impact upon the signal. 
The effect of these rejection criteria on the physics signal 
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FIG. 6: Distributions in P^'""'^ /pr, S, and Xr/E for show- 
ers from beam-halo muons (left) and photons (right) with 
Pt > 5 GeV/c. The beam- halo muon sample was obtained 
by requiring signals from the VW quadrant shadowing the 
reconstructed shower. The photon sample was obtained by 
imposing all the standard beam-halo muon rejection criteria 
with the exception of the requirement on the variable plotted. 



was checked using more restrictive selection criteria. The 
fraction of signal lost by the application of each of the 
muon-rejection criteria provided a measure of the sig- 
nal efficiency of these criteria. The application of the 
muon rejection criteria is estimated to result in the loss 
of f«7% of the direct photons at = 4 GcV/c, and 
Ri4% at pt = 7 GeV/c, for the 530 and 800 GeV/c beam 
data. The corresponding losses for the 515 GeV/c tt" 
beam data arc estimated to be 12% and 8% at = 4 
and 7 GeV/c, respectively. Corrections for these losses, 
which are functions of pt and t/cm, have been included in 
the results reported in this paper. 



E. Direct-photon candidates 

The most significant source of direct-photon back- 
ground is the electromagnetic decay of hadrons such as 
7r°'s and rj's. Photons were rejected from the direct- 
photon candidate sample if they combined with another 
photon in the same EMLAC octant to form a 77 pair with 
invariant mass, Mj^, in the 7r° or rj mass region and with 
energy asymmetry, A^^ = \E^^ — E^^\/{E^^ + Ej^), less 
than the value specified for each direct-photon candidate 
definition (described below). The residual background 
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FIG. 7: 77 mass distributions in the 7r° and rj (inset) mass 
regions for photon pairs with Aj^ < 0.75 and that satisfy the 
specified minimum px requirement in the data (histograms) 
compared to the results of the detailed detector simulation 
(•). The distributions have been normalized to the same area. 

from these mesons, as well as from other sources of back- 
ground, was evaluated using Monte Carlo samples that 
contained no generated direct photons. 

In the neutral meson analyses [10, 11], a 7r° candi- 
date was defined as a combination of two photons with 
A^^ < 0.75 detected within the fiducial region of the 
same EMLAC octant with 100 < M^^ < 180 MeV/c^; 
1] candidates were similarly defined with 450 < M^^ < 
650 MeV/c^ (Fig- 7). The energy asymmetry require- 
ment reduces uncertainties due to low energy photons. 
Background subtracted energy asymmetry distributions 
for photons in the mass region are shown in Fig. 8. 
Combinatorial backgrounds under the peak regions in the 
mass distributions were evaluated using sideband sub- 
traction. 

Three direct-photon candidate definitions were used 
in this analysis: 75A'', 90A'', and 755, which differed in 
the treatment of the background due to the decays of 
7r° mesons (the largest source of background). Correc- 
tions and backgrounds were calculated independently for 
each definition. In the 75iV (90iV) definition, any pho- 
ton that formed a combination with another photon in 
the same EMLAC octant with M^^ in the tt'^ mass re- 
gion and A~^^ < 0.75 {A^j < 0.9), was eliminated from 
the direct-photon candidate sample. A photon was also 
rejected if it formed a 77 pair with another photon in 
the same octant with invariant mass in the t] mass region 
with A^y < 0.8. 

The third definition, 755, rejected the same photons as 
the 75A^ definition. However, photons that formed mass 
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FIG. 8: Comparison of the side-band subtracted energy asym- 
metry distribution for photons from reconstructed 7r° decays 
in data (histogram) and the detailed detector simulation (•). 
Comparisons for two pr intervals, 4.0 < pr < 5.5 GeV/c and 
3.5 < pt < 4.0 GeV/c (inset) are shown. The distributions 
have been normalized to the same area. 



combinations with other photons in the tt*^ and rj side- 
band regions were included with a weight of 2 (when these 
photons were in the same octant) to explicitly account for 
direct-photon losses due to accidental combinations un- 
der the TT*^ and rj mass peaks. The 757V and 90N results 
both rely on the detailed simulation to account for the 
accidental losses of direct photons due to mass combina- 
tions with other photons that happen to land in the 
or 7] mass windows. 

Each of these definitions have relative strengths and 
weaknesses. The 90A^ definition rejects the most back- 
ground, however, the residual background in this defini- 
tion is more sensitive to the simulation of the low energy 
photons from tt" decay, which cause the losses at high 
Ajj. Additionally, accidental direct-photon losses are 
largest for this definition and must be accounted for by 
the simulation. The 75iV and 755 definitions are less 
sensitive to the details of the low energy photon simu- 
lation, but lead to larger background levels. Comparing 
the fully corrected results using three definitions provided 
insights into the systematic uncertainties associated with 
the direct-photon cross-section measurements. 



F. Detailed detector simulation 

The Meson West spectrometer was modeled using a de- 
tailed GEANT [28] simulation (DCS). The simulated par- 
ticles in the generated events were propagated through 



the GEANT model of the spectrometer and interacted ap- 
propriately. Simulated showers in the EMLAC were pa- 
rameterized once the energies of the interacting particles 
were <10 MeV. A preprocessor converted geant infor- 
mation into the hits and strip energies simulating the 
detailed response of the various detectors. The prepro- 
cessor simulated hardware effects such as channel noise 
and gain variations. 

Single particle distributions, reconstructed data 
events, and the herwig v5.6 [29] event generator were 
employed as inptits to the detailed geant simulation. 
The direct-photon backgrotmd studies were carried out 
primarily using samples of HERWIG 2—^2 QCD hard 
parton scatters (IPR0C=1500). No direct photons were 
included in these samples. The detailed simulations were 
used to evaluate contributions to the direct-photon back- 
ground from numerous sources including 7r°, rj, lu, tj', K1, 
Kg, e^, neutrons and charged pions. 

Generated events were processed through the recon- 
struction software used for the data analysis. This tech- 
nique accounted for inefficiencies and biases in the re- 
construction algorithms. Additional information on the 
detailed simulation of the Meson West spectrometer can 
be found elsewhere [10, 30]. 

The simtilated herwig events were weighted (in an it- 
erative fashion) so that the resulting tt*^ spectra agreed 
with our measured cross sections as functions of px and 
Ucm [10, 11]. As a result, the final corrections were based 
on the shape of the spectra from data rather than on the 
HERWIG event generator [30]. Similarly, the generated 
rj's were weighted so that the resulting simulated produc- 
tion spectra were consistent with our measured ri/w^ re- 
sults [31]. HERWIG v5.6 was assumed to reproduce other 
particle spectra after these were weighted by the results of 
the fits to the 7r° data [32, 33] . Photons reconstructed in 
the background samples that satisfied the direct-photon 
candidate requirements, ^bckg, were assigned the weight- 
ing appropriate to their identified parent. 

Figure 7 shows 77 mass spectra in the and rj mass 
regions compared to the results of the detailed GEANT 
simulation. In addition to demonstrating that the DGS 
simulated the EMLAC resolution well, the agreement be- 
tween the levels of combinatorial background indicates 
the DGS provided a reasonable simulation of the un- 
derlying event structure. Figure 8 shows a comparison 
between the DGS results and the data for the 7r° energy 
asymmetry after side-band subtraction. This figure illus- 
trates that the simulation accurately describes the losses 
of low energy photons. 

The direct-photon reconstruction efficiency was eval- 
uated using a dedicated sample of HERWiG-generated 
direct-photon events. Reconstuction inefficiencies for di- 
rect photons were relatively small over most of the kine- 
matic range (Fig. 9). 
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FIG. 9: Probabilities to reconstruct direct photons using the 
755, 75A'^, and 90N definitions as functions of generated pr 
of the photons for several rapidity ranges as determined via 
the detailed GEANT simulation. 
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FIG. 10: Ratios of reconstructed to generated energy for tt^'s 
and direct photons as functions of their generated pr values 
as determined via detailed GEANT simulations. The small 
offset from 1.0 is due to the combined effects of the steeply 
falling particle production spectra and the EMLAC resolu- 
tion. 



eratcd energy as functions of generated px for tt'^'s and 
direct photons from the DGS. Cuts were applied to both 
the photons from 7r° decay and direct photons ensur- 
ing equal treatment by the reconstruction software [34]. 
These two distributions agree well. 



H. Direct-photon background 



G. EMLAC energy scale 

The calibration of the EMLAC response was based 

on the reconstructed masses of 7r° mesons in the 77 de- 
cay mode and verified using yy's, lo's, converted photons 
(e+e"), and J/ip's [19, 34]. The steeply falling 7r° pr 
spectrum, combined with the calorimeter's resolution, re- 
sults in a small offset (1%) in the mean reconstructed 
photon energies when plotted as a function of the gener- 
ated pt (Fig- 10). To account for this offset and for po- 
tential biases in the calibration procedure, the simulated 
EMLAC was calibrated in the same manner as the real 
detector. The impact of detector resolution on the en- 
ergy scale calibration and on the production spectra was 
incorporated into the Monte Carlo based corrections [30] . 
The simulation was also employed to evaluate the mean 
correction for energy deposited in the material upstream 
of the EMLAC. Two-photon mass distributions in the 7r° 
and T] mass regions from both data and DGS events after 
energy scale calibration are illustrated in Fig. 7. 

The energy scale calibration was also verified to be ap- 
propriate for direct photons. This is important since the 
EMLAC energy response was sensitive to event struc- 
ture; direct photons are expected to be more isolated 
than corresponding photons from meson decays since the 
mesons arc frequently accompanied by other jet particles. 
The single photon energy scale was verified with 77 ^ 77 
where the photons were widely separated, as well as with 
isolated electrons, and the DGS. Figure 10 displays a 
comparison of the ratio between reconstructed and gen- 



The residual background to the direct-photon signal 
from tt'^'s and r^'s, as well as from other (previously men- 
tioned) sources, was calculated by applying the 75A^, 
90iV, and 755 photon definitions to DGS samples that 
contained no generated direct photons. Ratios of the re- 
sulting ^bckg to the measured tt*^ cross sections [10, 11] are 
shown as functions of pr for interactions on the beryllium 
target in Fig. 11. As indicated in the section describing 
the various direct-photon candidate definitions, the 90N 
background level is substantially lower than the 75A'' and 
7519 background levels. 

Figure 12 shows the jtckg/'^'^ for several ycm inter- 
vals for the simulated 530 GeV/c pBe interactions. The 
shapes of the background are quite similar in the back- 
ward and central ycm regions. However, at forward ycm 
and high pT the contribution to the background due to 
coalesced photons from tt'^ decays is clearly visible. The 
Jbckg/i^'^ were fitted to functions in pT and j/cm- Projec- 
tions from the fits are indicated by the curves in Figs. 11 
and 12. 



I. Parameterized detector simulation 

A second Monte Carlo simulation of detector effects 

(PMC) was used to cross check the DGS and for large 
statistics studies. This simulation employed parame- 
terizations of physics cross sections and detector re- 
sponses [30, 34]. The inclusive tt" and direct-photon cross 
sections were parameterized as two-dimensional functions 
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FIG. 11: 76cfcg/7r'' for each direct-photon candidate definition 
as a function of pr based upon analysis of DGS samples ex- 
cluding the direct-photon samples. The points represent the 
binned results; the curves represent fits to these backgrounds 
integrated over rapidity. 



FIG. 13: Top: comparison between background-subtracted 
77 meiss distributions for data (•) and the PMC (histogram) 
for photon pairs with A-,^ < 0.75 in the n'^ and r; mass re- 
gions. Note that the background levels are centered on zero 
(dashed line) in the data mass plots due to the background 
subtraction. Bottom: comparison of the corresponding en- 
ergy asymmetry distributions for 7's from tt" decays. The 
distributions have been normalized to the same area. 
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FIG. 12: jbckg/T^^ for each direct-photon candidate defini- 
tion as a function of pr for backward, central, and forward 
rapidities based upon analysis of DGS samples excluding the 
direct-photon samples. The curves represent the results of the 
background fits integrated over the indicated rapidity range. 



ill px cind ycm 

[30]. The T], u!, and 77' cross sections were 
represented using values consistent with our measured 
7?/7rO [10, 11, 31], w/ttO [13], and using ?777rO=0.85, as 
indicated in Ref. [35] (after accounting for the updated 
1]' 77 branching ratio [36]). 

An interaction vertex was generated in the simulated 
target for every event (Fig. 3). This primary vertex dis- 
tribution was used to ensure the simulated photons tra- 
versed the proper amount of material. Generated mesons 
were allowed to decay into final state particles. Photons 
were allowed to convert into e+e~ pairs when appropri- 
ate; the energy of the resulting electrons was reduced 
using the GEANT function for bremsstrahlung radiation. 
Electron four-vectors were smeared for multiple scatter- 
ing in the target region and the resolution of the tracking 
system, and adjusted to account for the impulse due to 
the magnet. The parameters of the photons that did 
not convert were smeared to account for energy and po- 
sition resolution [19]. Figure 13 displays a comparison 
between 77 mass distributions from the PMC and the 
background-subtracted data in the 7r° and rj mass regions 
as well as the energy asymmetry distribution for photons 
from 7r° decays. As illustrated by these plots, the PMC 
provides a reasonable characterization of the data. 

The PMC was used to identify contributions to the 

direct-photon background. Photons from tt^ > 77 

decays with A^^ greater than required by the relevant 
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FIG. 14: Contribution to ^bckg/T^^ from tt" decays for the 
755 photon definition evaluated via the PMC. Also shown 
aie the contributions from each of the categories described in 
the text. 

photon definition were the dominant contribution to the 
background. However, photons from tt" decays with 
smaller A^^ also contribute to ^bckg when the other pho- 
ton is lost. Hence, the remainder of the photons from 7r° 
decays were classified according to otir ability to detect 
both photons from the 7r° decay. The photons were pro- 
jected to the front face of the EMLAC. If one was outside 
the EMLAC fiducial volume while the other was within 
that fiducial volume, then the latter contributed to jtckg- 
In addition, photons that landed in different EMLAC oc- 
tants were classified as background candidates. Contri- 
butions to jbckg from photon conversions and detector 
inefficiencies were obtained from the sample when both 
photons landed within a single EMLAC octant. Each 
photon was assigned a weight based on the conversion 
probability of the other photon. The 7t,cfeg contribution 
from detector inefficiencies was determined by weight- 
ing each photon by the non-detection probability for the 
other photon. Double counting between these sources 
was eliminated by applying the non-conversion probabil- 
ity of the undetected photon to each jbckg contribution. 
Since the "fbckg were required to be inside the fiducial vol- 
ume, the impact of the EMLAC acceptance was properly 
taken into account. 

The DGS was used to assess the coalescence proba- 
bility function as the coalescence of two photons from a 
tt'^ decay into a single reconstructed shower depended on 
many factors including the relative geometry of the two 
photons, their total and relative energies, and the ability 
of the reconstructor to resolve the photons. Background 
contributions from this source were assigned the 7r° en- 



FIG. 15: Contribution to 76cfcg/7r*' from 7r°'s, ?7's, w's, and r7''s 
for the 75S photon definition as evaluated using the PMC. 



ergy and weighted by the non-conversion probability and 
detection efficiency of each photon and the acceptance of 
the higher energy photon. 

The ratio of ^bckg from tt" 77 decays to generated 
7r°'s is shown in Fig. 14 for incident 530 GeV/c p beam 
on the beryllium targets. The contributions from each of 
the categories described above are also shown in the same 
figure. It is evident from the figure that photons result- 
ing from highly asymmetric 7r° decays are the dominant 
contribution to the direct-photon background, although, 
at very high px, the coalescence contribution is also sig- 
nificant. 

The background contribution from the decay 77 — > 77 
was evaluated in an analogous manner. Photons from 
r] decays with energy asymmetry Aj^ > 0.8 were in- 
cluded in the direct-photon background. Photons from 
the decays to — > vr''7, 77' p^, rj' uj^f, rj' 77, and 
7r° 7e"'"e~ were also evaluated as potential contribu- 
tors to 76cfcg- Photons arising from the subsequent decay 
of the 77° in the case of the u decay and the p — > 7r°7r° de- 
cays were not considered since those contributions were 
already included in the simulation of the inclusive 7r° 
sample. The background contribution from the e+e~ 
pair in 7e+e~ decays was negligible due to the 

small branching ratio for this decay mode combined with 
the typically significantly lower px of the daughter elec- 
trons relative to the parent 7r°. Furthermore, most elec- 
trons were rejected by the distance-to-nearest-track re- 
quirement (Fig. 5). 

The total 'jbckg/'^^ as determined using the PMC is 
shown in Fig. 15 for incident 530 GeV/c p on beryllium. 
The contribution to ^bckg/T^'^ due to each of the particles 
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FIG. 16: Diflcrcncc between ^bckg/i^'^ for the hydrogen and 
beryUium targets as a function of pr as determined via the 
PMC. The dashed curve represents a fit to this difference. 

considered in the PMC is also shown in that figure. The 
contribution from 77 decays is roughly 20% of the tt'^ con- 
tribution, as expected from simple considerations of the 
T] and tt" production rates [10, 11] and their respective 
two-photon branching ratios [36]. 

J. Direct-photon background for hydrogen target 
data 

Fits to ^bckg/T^^ evaluated using the DGS were only 
made for the beryllium target due to relatively poor DGS 
statistics in the other targets. However, ^bckg/T^^ is ex- 
pected to be slightly different for each target due to the 
different amounts of target material the photons must 
traverse. A correction to ^hckg/T^^ was therefore eval- 
uated using the PMC. The difference between ^bckg/i^^ 
on the hydrogen and beryllium targets is shown for the 
800 GeV/c proton beam sample in Fig. 16. The corre- 
sponding results are similar for the other incident beams. 
Note that the difference is small (< 0.003) but positive at 
low pt because the conversion background contribution is 
larger for the more upstream hydrogen target. However, 
at high Pt, the correction becomes very slightly negative 
because the contribution due to coalescence is reduced 
by the additional photon conversions. The target differ- 
ences were fit as functions of pT and ycm for each incident 
beam and applied as additive corrections to the nominal 
IhckglT^^ fit obtained using the DGS and the beryllium 
target. 

K. Background subtraction uncertainties 

The systematic uncertainty in the background subtrac- 
tion was estimated via the sensitivity of the results to 
the various direct-photon definitions and to variations in 
the size of the ix^ and ry side-band regions and via de- 
tailed studies of the Monte Carlo simulations. Ratios 
of direct-photon cross sections obtained using the differ- 
ent candidate definitions arc displayed in Fig. 17 for all 
three incident-beam samples. These ratios generally dif- 
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FIG. 17: Ratios between direct-photon cross sections evalu- 
ated using the 90N and 75S direct-photon candidate defini- 
tions as functions of pliotoii pr for the three major incident 
beams. Ratios between direct-photon cross sections evaluated 
using the 75N and 75S direct-photon candidate definitions are 
also shown. 

fered from unity by < 5%. This very good agreement 
between the direct-photon cross sections measured using 
these different definitions illustrates the robust nature of 
our results. 

Comparisons between detailed and parameterized de- 
tector simulations provided increased confidence in our 
understanding of the direct-photon background since the 
two simulations represent essentially independent calcu- 
lations of jbckg ■ The PMC-evaluated -ybckg is based on rel- 
atively simple considerations and included only the major 
background sources. Special runs of the DGS were per- 
formed limiting the background sources to those included 
in the PMC. "fbckg/'^'^ from the restricted background 
source DGS and the PMC are compared in Fig. 18. The 
two simulation results agree quite well. 

Measured four-vectors from particles reconstructed in 
data events were also employed as inputs to the GEANT 
simulation. This provided an additional measure of the 
uncertainty associated with the modeling of the detec- 
tor environment. A comparison of Jbckg/T^^ produced 
in the DGS using the weighted herwig event generator 
and reconstructed data events is shown in Fig. 19. Only 
background sources common to both simulations were in- 
cluded in the comparison. The jbckg/T^^ results evaluated 
via these two methods are consistent. 

The impact of uncertainty in 77/71° on jbckg/T^'^ was 
determined using the PMC. The values used for rj/ir^ 
were shifted by ±0.04; these shifts resulted in a ±1.5% 
change in jbckg/T^°- 
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FIG. 19: Comparison of jbckg/'!T° from the DGS using 
weighted herwig events and using reconstructed data inputs 

as the event generators for the 755, 75N, and 90N direct- 
photon candidate definitions. 



L. Normalization 

Electronic scalers that counted signals from the beam 
hodoscope and beam hole counters were used to estab- 
lish the number of beam particles incident on the target. 
Other scalers logged the state of the trigger and data 
acquisition system. Information from these scalers was 
used to determine the number of beam particles incident 
upon the spectrometer when it was ready to record data. 
Small corrections for multiple occupancy in the beam ho- 
doscope and for absorption of beam in the target material 
were included in the normalization. 

Additional discussion on the normalization includ- 
ing details regarding independent normalization methods 
and cross checks can be found elsewhere [10, 11]. The 
net systematic uncertainty in the overall normalization 
is 8%. This evaluation is based upon several considera- 
tions: the good agreement between results from the inde- 
pendent normalization methods, the stability of the cross 
section results from various sections of the data, the in- 
ternal consistency of the scalers, and a detailed analysis 
of the design, implementation, and performance of the 
trigger. 



M. Secondary beam contamination 

The Meson West beamline transported 0.8 TeV/c pri- 
mary protons from the Fermilab Tevatron and 0.5 TeV/c 
secondary beams of both polarities. The 800 GeV/c pri- 



mary beam has a negligible momentum spread. The 
negative secondary pion beam was determined to have 
a mean momentum of 515 ± 2 GeV/c with an estimated 
halfwidth of 30 GeV/c. The secondary proton beam was 
determined to have a mean momentum of 530 ± 2 GeV/c 
and similar momentum spread. This momentum spread 
in the secondary beams introduces a small uncertainty 
(< 4%) in comparisons with theory evaluated at the 
mean beam momentum. 

The minority particle fractions in the 530 GcV/c pos- 
itive secondary beam were measured using the beamline 
Cerenkov detector to be 2.75% n+ and 0.5% K+ [24]. 
Analysis of the 515 GeV/c negative secondary beam in- 
dicated 1% K~ component, and negligible p content. 

The direct-photon cross sections for incident 
530 GeV/c p beam have been corrected for the small 
incident pion contamination as follows. An estimate of 
the direct-photon cross section from incident 7r+ was 
generated from a product of our measured incident 7r~ 
cross sections and the ratio of next-to-leading-order 
PQCD calculations [37, 38] for direct-photon production 
in TT"*" and tt~ beams to account for the differences 
in incident quark content. The direct-photon cross 
sections reported here for incident 530 GeV/c protons 
were corrected to account for the minority pions using 
this estimate of the direct-photon cross section from 
incident tt"*" combined with the measured incident 
minority particle fraction. The correction was < 1% at 
Pt = 4 GeV/c and < 5% at pt = 8 GeV/c. 

The direct-photon cross-section measurements have 
not been corrected for the incident kaon contributions 
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FIG. 20: Relative systematic uncertainty for the direct- 
photon cross section as a function of pr for incident 
530 GeV/c protons. 



due to the lack of an incident kaon PDF. The potential 
bias from this source was studied using restrictive cuts 
based upon the Cerenkov detector data and the impact of 
the kaon contamination upon the reported cross-section 
measurements is expected to be negligible. 



N. Summary of systematic uncertainties 

The principal contributions to the systematic uncer- 
tainty in the direct-photon cross-section results arose 
from the following sources: normalization, background 
subtraction, the calibration of photon energy response 
and detector-resolution unsmearing, and incident beam 
contamination (for the 530 GeV/c p beam). Other 
sources of uncertainty, which contributed at the one 
to two percent level, included: photon reconstruction, 
beam-halo muon rejection, geometric acceptance, photon 
conversions, trigger response, and vertex finding. The 
relative systematic uncertainty in the direct-photon cross 
section is presented as a function of pt in Fig. 20 for 
the 530 GeV/c pBe data. With the exception of the 
beam contamination contribution, the corresponding un- 
certainties for the other measurements are similar. The 
total systematic uncertainties, combined in quadrature, 
are quoted with the cross sections in the appropriate ta- 
bles. Note that some of these contributions to the sys- 
tematic uncertainty (e.g. normalization) arc strongly cor- 
related between bins. 



FIG. 21: Ratios of the 90N direct-photon candidate spectra to 
the measured tt*^ cross sections (points) compared to ■jbckg/T^'^ 
from the DGS (curves) as functions of pt for the data samples 
considered in this paper. The error bars represent statistical 
contributions to the uncertainties. The width of each band 
around the background represents the systematic uncertainty 
on that background. 
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FIG. 22: Invariant differential cross sections per nucleon for 
direct-photon production as functions of pr, averaged over 
rapidity, for 515 GeV/c n~ and 800 and 530 GeV/c proton 
beams incident upon beryllium. The error bars represent the 
statistical and systematic uncertainties combined in quadra- 
ture; the innermost interval indicates the statistical uncer- 
tainties. 
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FIG. 23: Ratios of inclusive direct-photon production per nu- 
cleon in pBe to those in pp collisions as functions of photon 
Pt, compared with curves representing predictions from Hl- 
JING. The error bars represent only statistical contributions 
to the measurement uncertainties. 
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FIG. 24: Invariant differential cross sections per nucleon for 
direct-photon production as functions of pr in 7r~Be and iv^p 
interactions at 515 GeV/c. The error bars represent the sta- 
tistical and systematic uncertainties combined in quadrature; 
the innermost interval indicates the statistical uncertainties. 
Overlaid on the data are NLO PQCD results for three choices 
of the factorization and renormalization scales. 



FIG. 25: Invariant differential cross sections per nucleon for 
direct-photon production as functions of pr in pp interactions 
at 800 and 530 GeV/c. The error bars represent the statistical 
and systematic uncertainties combined in quadrature; the in- 
nermost interval indicates the statistical uncertainties. Over- 
laid on the data are NLO PQCD predictions for the CTEQ5M 
and MRST2003C PDF; the band illustrates the PDF uncer- 
tainty estimated via the MRST2001E set. 



IV. RESULTS AND DISCUSSION 
A. Results 

The ratios of the 90A'' direct-photon candidate spectra 
to the measured tt*^ cross sections [10, 11] as functions of 
Pt are compared to jbckg/'^''^ from the DGS in Fig. 21. All 
three samples are characterized by rapidly rising direct- 
photon signal to background ratios. 

The direct-photon cross sections were extracted statis- 
tically using the fits to ^bckg/T^^ ■ The invariant differen- 
tial cross sections per nucleon for direct-photon produc- 
tion from 515 GeV/c tt^ and 800 and 530 GeV/c pro- 
ton beams incident on beryllium are presented as func- 
tions of Pt in Fig. 22. Results from 515 GeV/c 7r~ and 
530 GcV/c p beams are averaged over the rapidity range 
—0.75 < j/cm < 0.75; results from the 800 GeV/c p beam 
are averaged over —1.0 < j/cm < 0.5. Data points are 
plotted at abscissas that correspond to the average value 
of the cross section in each pT bin, assuming local ex- 
ponential dependence on pT [39]. The inclusive cross 
sections are tabulated in the Appendix (Tables I — VIII) . 

Ratios of the direct-photon cross sections per nucleon 
on beryllium to those on hydrogen are shown versus pt 
in Fig. 23. The ratios are compared with results from 
the Monte Carlo program hijing [40]. The results from 
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FIG. 26: The ratio of 800 GeV/c to 530 GeV/c proton beam 
direct-photon cross sections as a function of xt- The er- 
ror bars represent the statistical and systematic uncertain- 
ties combined in quadrature; the innermost interval indicates 
the statistical uncertainties. The dotted curve shows the 
resuhs of LO PQCD calculations using the CTEQ5L PDF 
and fi — Pt/2. Shaded bands illustrate the results from 
NLO PQCD calculations. The outer band is bounded by 
M = Pt/3 and jj, — pr; the inner band by the uncertainty 
of the MRST2001E PDF set. 



HIJING, which is designed to simulate particle production 
in pp, pA, and AA collisions, are in good agreement with 
the data. The PQCD calculations described below for 
interactions on beryllium have been adjusted to account 
for nuclear effects using hijing results. 



B. Comparisons with PQCD 

Next-to-leading order (NLO) PQCD calculations [37, 
38] using the GRV92 PDF [7] are compared to our mea- 
sured direct-photon cross sections for incident 7r~ beam 
in Fig. 24. Theoretical results are presented for three 
choices of factorization scale: /i — pt, VtI^^ and pr/'i- 
The renormalization scale was set equal to the factoriza- 
tion scale. Expectations for these choices of /i lie signif- 
icantly below the data, in addition to exhibiting a large 
dependence on choice of scale. Calculations that include 
higher-order effects, such as threshold resummation and 
certain NNLO diagrams [41-45] significantly reduce scale 
dependence relative to NLO calculations, and yield re- 
sults comparable to the NLO calculation with = pt/2. 

There is additional theoretical uncertainty associated 
with the choice of PDF. Recently published PDF pro- 
vide individual fits for variations in fit parameters. Un- 



FIG. 27: The ratio of 515 GeV/c tt" beam to 800 GeV/c 
proton-beam direct-photon cross sections as a function of pr- 
The error bars represent the statistical and systematic un- 
certainties combined in quadrature; the innermost interval 
indicates the statistical uncertainties. Results of LO PQCD 
calculations using /i = pt/2 are represented by the dotted 
curve. The shaded band illustrates results from NLO PQCD 
calculations bounded by /i = pt/3 and fi = pr- 



certainty bands based upon these fits are presented in 
Fig. 25 for the MRST2001E PDF set [46, 47] compared 
with the pp data. We conclude there is insufficient flexi- 
bility in the PDF to accommodate the difference between 
theory and data. 

Not surprisingly, NLO PQCD has reduced scale depen- 
dence and significantly improved agreement with data 
when compared to ratios of direct-photon cross sections. 
The ratio of the 800 to 530 GeV/c direct-photon cross 
sections is presented as a function of xt = '2pt/^/s in 
Fig. 26. The ratio between 515 GeV/c tt~ beam and 
800 GeV/c proton-beam direct-photon cross sections as 
a function of is shown in Fig. 27. The broad bands on 
the NLO theory in both figures are bounded by /i — ptI^ 
and pt\ the inner band in Fig. 26 represents the uncer- 
tainty attributed to the PDF. The dotted curves in both 
figures show the results of leading-order (LO) PQCD cal- 
culations for /i — pt/'2,. We chose xt as the plotting 
variable in Fig. 26 to compensate for the different aver- 
age parton-parton collision energies in the numerator and 
the denominator; pT was used as the plotting variable 
in Fig. 27 since the average energy per colliding valence 
quark is already similar in the numerator and the denomi- 
nator. There is relatively little difference between the LO 
and NLO curves in Fig. 26, where the beam particle types 
are the same for both the numerator and the denomina- 
tor. In marked contrast, however, in Fig. 27, where the 
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FIG. 28: Invariant differential cross sections per nucleon for 
direct-photon production as functions of pr, averaged over 
rapidity, for 515 GcV/c 7r~ and 800 and 530 GcV/c proton 
beams incident upon beryllium. The error bars represent the 
statistical and systematic uncertainties combined in quadra- 
ture; the innermost interval indicates the statistical uncertain- 
ties. Overlaid on the data are NLO PQCD and fcr-enhanced 
NLO PQCD calc:ulations. GRV92 PDF were used in the inci- 
dent calculations, while CTEQ5M PDF were used in the 
incident proton calculations. 



beam particle types arc different in the numerator and 
the denominator, the LO and NLO calculations are very 
different, and only the NLO result is consistent with our 
data for the PDF employed in these calculations. 

Our direct-photon cross-section measurements, as il- 
lustrated in Figs. 24 and 25, as well as our 7r° cross- 
section measurements reported in previous publica- 
tions [9-11], are not satisfactorily described by the avail- 
able NLO PQCD calculations. Similar discrepancies have 
been observed between PQCD calculations and other 
measurements of high-pT direct-photon and tt" produc- 
tion [48-50] and are attributed to the effect of initial- 
state soft-gluon radiation (sec also Rcfs. [51-53]). The 
importance of including gluon emission through the re- 
summation formalism has long been recognized, and the 
corresponding calculations have been available for some 
time for the diphoton process [54, 55]. However, these 
calculations have only recently been considered for inclu- 
sive direct-photon production [56-59]. While these de- 
velopments are encouraging (e.g., see Fig. 2 of Ref. [58]), 
the calculations are still incomplete. In the absence of 
a rigorous treatment, the phenomenological prescription 
described in Refs. [9, 48] is used in the following com- 
parisons of PQCD with our data. In this approach, 
soft-gluon radiation is parameterized in terms of an ef- 



FIG. 29: Invariant differential cross sections per nucleon for 
direct-photon production as functions of j/cm in pBe interac- 
tions at 530 GcV/f for several intervals in pT- The error bars 
represent the statistical and systematic uncertainties com- 
bined in quadrature; the innermost interval indicates the sta- 
tistical uncertainties. Overlaid on the data are fcT-enhanced 
NLO PQCD calculations using the CTEQ5M PDF. 



fective (kr) that provides an additional transverse im- 
pulse to the incoming partons. Since the hard-scattering 
cross sections fall steeply as functions of pt, the net ef- 
fect of this additional transverse impulse is to enhance 
high-pT yields. The kr enhancement factors (as func- 
tions of Pt) are evaluated by computing ratios of results 
from LO PQCD calcidations [3] with and without fcr- 
These fcy enhancement factors are applied to the results 
of the NLO PQCD calculations. The (kx) values used 
in generating the kx enhancement factors arc represen- 
tative of values that describe kinematic distributions ob- 
served in production of high-mass tt'^tt'^, 'Jtt^, and 77 
systems [9, 34]. The fcr-cnhanced NLO calculations are 
compared to our direct-photon measurements in Fig. 28, 
and provide a much improved description of our data 
(relative to the NLO calculations). Comparisons of the 
/cT-enhanced NLO calculations to our data as functions 
of 2/cm for several pr intervals are shown in Fig. 29 for 
530 GeV/c p beam on berylHum. 

Results of kr enhanced NLO PQCD calculations using 
CTEQ6.1M PDF [60, 61] are compared to our direct- 
photon measurements for incident protons in Fig. 30. 
The CTEQ6.1 PDF uncertainty bands are also shown in 
that figure. While the uncertainties associated with the 
PDF are small at low Xt, they are significantly larger 
than the uncertainties in data at moderate and high xt, 
and are dominated by the uncertainty in the gluon distri- 
bution. There is also a large variation in the central value 
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FIG. 30: Invariant differential cross sections per nucleon for 
direct-photon production as functions of xt in pBe interac- 
tions at 800 and 530 GeV/c. The error bars represent the 
statistical and systematic uncertainties combined in quadra- 
ture; the innermost interval indicates the statistical uncer- 
tainties. Overlaid on the data are fcT-enhanced NLO PQCD 
calculations. The shaded region represents the uncertainty 
band associated with the CTEQ6.1M PDF set. 



of several PDF sets, as illustrated in Fig. 31. The ma- 
jor difference between the CTEQ5M [62] and CTEQ6.1M 
PDF is that CTEQ6.1M has a much harder gluon (moti- 
vated by measurements of the inclusive jet cross section 
at the Tevatron collider). Cross sections calculated us- 
ing MRST2003C PDF [63] are similar to results using 
CTEQ5M PDF. Figure 32 shows the ratio of the direct- 
photon cross sections from 800 GeV/c and 530 GeV/c 
p beams on beryllium, as a function of pT^ compared 
to results from A:T-enhanced NLO calculations using the 
MRST2003C, CTEQ5M, and CTEQ6.1M PDF. Calcula- 
tions with PDF that have softer gluons provide a better 
description of our data at high pT- 



V. CONCLUSIONS 

High-pT direct-photon production has been measured 
in interactions of 515 GeV/c tt" and 530 GeV/c and 
800 GeV/c protons with beryllium and hydrogen targets. 
The inclusive direct-photon cross sections were compared 
with NLO PQCD and fcT-enhanced NLO PQCD calcu- 
lations for several choices of PDF. Without kr enhance- 
ment, available NLO PQCD calculations do not satisfac- 
torily represent our data. The data are described better 
by kx enhanced NLO PQCD calculations using the softer 
gluons of CTEQ5 and MRST2003 than the harder gluons 



FIG. 31: Invariant differential cross sections per nucleon for 
direct-photon production as functions of xt in pBe interac- 
tions at 800 and 530 GeV/c. The error bars represent the 
statistical and systematic uncertainties combined in quadra- 
ture; the innermost interval indicates the statistical uncer- 
tainties. Overlaid on the data are fcT-enhanced NLO PQCD 
calculations using CTEQ5M and CTEQ6.1M PDF. 
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FIG. 32; The ratio of 800 to 530 GeV/c proton-beam direct- 
photon cross sections as a function of pr- The error bars rep- 
resent the statistical and systematic uncertainties combined 
in quadrature; the innermost interval indicates the statistical 
uncertainties. Results from fcr-enhanced NLO PQCD calcu- 
lations are overlaid on the data. 
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of CTEQ6.1. As mentioned in the Introduction, direct- 
photon measurements provide potential constraints on 
the gluon, particularly at large x. This is especially rele- 
vant since recent global fits have used jet data from Teva- 
tron Run I to constrain the large-x gluon — thereby pos- 
sibly incorporating new physics effects into the PDF [61]. 
We note that recent studies have suggested potential in- 
compatibilities between the small-./; gluon determined us- 
ing HERA data and the large- .t gluon required to fit the 
Tevatron jet data [64]. Inclusion of direct-photon data 
in the global fits would provide additional constraints on 
the gluon distribution, independent of the jet data from 
the Tevatron, thereby enhancing the discovery potential 
of the jet data being acquired in Tevatron Run II. 
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APPENDIX: TABULATED CROSS SECTIONS 



TABLE L Invariant differential cross sections {Eda/d^p) per nucleon for direct-photon production in pBe collisions at 800 and 
530 GeV/c, and 7r~Be collisions at 515 GeV/c as functions of pt- 



Pt 


prJe at oUU /c 


prSe at boU OeV/c 
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— l.U S Vcm S U.O 
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— U. /O S 2/cm S U./O 


(GeV/c) 


|nb/(GeV/c) J 


[nb/^OeV/cj J 


[lib / [KjCV / c) \ 


3.50 - 


3.75 


2.99 ± 0.24 ± 0.62 


1.83 ± 0.17 ± 0.37 


1.79 ± 0.10 ± 0.33 


3.75 - 


4.00 


■1 f f\f\ 1 n r\f\ "1 1 n orv 

1.500 ± 0.091 ± 0.29 


0.993 ± 0.062 ± 0.18 


0.924 ± 0.036 ± 0.16 


4.00 - 


4.25 


0.908 ± 0.031 ± 0.16 


0.478 ± 0.019 ± 0.081 


0.5058 ± 0.0081 ± 0.080 


4.25 — 


4.50 


r\ A an 1 r\ f\-t t i f~\ r\>'7f\ 

0.482 ± 0.011 ± 0.079 


0.2508 ± 0.0047 ± 0.040 


0.2891 ± 0.0052 ± 0.043 


4.50 


4.75 


0.2679 ± 0.0068 ± 0.041 


0.1334 ± 0.0030 ± 0.020 


0.1736 ± 0.0034 ± 0.024 


4.75 
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0.1616 ± 0.0042 ± 0.023 


0.0757 ± 0.0021 ± 0.011 


0.1084 ± 0.0024 ± 0.014 
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5.00 - 


5.25 


98.4 ± 3.1 ±14 


41.3 ± 1.4 ±5.5 


70.8 ± 1.8 ±9.0 
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63.7 ± 2.3 ±8.4 


24.4 ± 1.0 ±3.2 


42.7 ± 1.3 ±5.2 
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20.02 ± 0.77 ± 2.4 
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12.86 ± 0.52 ± 1.6 


3.98 ± 0.25 ±0.49 


10.38 ± 0.37 ± 1.2 


6.50 - 
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5.36 ± 0.30 ±0.63 


1.79 ± 0.15 ±0.22 


4.74 ± 0.23 ±0.55 
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7.50 


2.51 ± 0.19 ±0.29 


0.718 ± 0.086 ± 0.089 


2.12 ± 0.15 ±0.25 
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1.11 ± 0.12 ±0.13 


0.324 ± 0.054 ±0.041 


1.148 ± 0.099 ±0.14 
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0.420 ± 0.052 ± 0.049 
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10.00 - 
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0.00095 ± 0.00095 ± 0.00015 


0.0060 ± 0.0028 ± 0.0008 



TABLE II: Invariant differential cross sections [Eda/d^p) for direct-photon production in pp collisions at 800 and 530 GeV/c, 
and ■K~p collisions at 515 GeV/c as functions of pr- 

p^ pp at 800 GeV/c pp at 530 GeV/c tt'p at 515 GeV/c 

-1.0 < J/cm < 0.5 -0.75 < J/cm < 0.75 -0.75 < 7/cm < 0.75 

(GcV/c) [nb/(GeV/c)^] [nb/(GcV/c)^] [nb/(GeV/c)^] 

3.50 - 4.00 1.88 ± 0.30 ±0.38 1.16 ± 0.19 ± 0.22 1.72 ± 0.42 ± 0.31 

4.00 - 4.50 0.600 ± 0.044 ±0.10 0.375 ± 0.030 ± 0.062 0.493 ± 0.028 ± 0.076 
4.50 - 5.00 0.2255 ± 0.0094 ±0.034 0.1076 ± 0.0045 ±0.016 0.165 ± 0.011 ± 0.022 

[pb/(GeV/c)-] [pb/(GeV/c)-] [pb/(GeV/c)^^ 

5.00 - 5.50 82.2 ±4.4 ±11 33.6 ±2.2 ±4.4 64.8 ±6.0 ±8.1 

5.50 - 6.00 35.0 ±2.4 ±4.4 11.4±1.1± 1.4 25.5 ±3.5 ±3.0 

6.00 7.00 9.80 ± 0.74 ± 1.2 3.43 ± 0.37 ± 0.42 9.0 ±1.3 ±1.0 

7.00 - 8.00 1.73 ± 0.28 ±0.20 0.42 ± 0.11 ± 0.05 2.81 ± 0.62 ± 0.33 

8.00 - 10.00 0.339 ± 0.070 ± 0.039 0.010 ± 0.015 ± 0.001 0.30 ± 0.14 ± 0.04 
10.00 - 12.00 0.017 ± 0.014 ± 0.002 
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TABLE III: The invariant differential cross section {Eda/dPp) per nucleon as a function of j/cm and pr for the inclusive reaction 
pBe 7X at 530 GeV/c. 



Pt -0.75 < 7/cm < -0.50 -0.50 < 7/cm < -0.25 -0.25 < i/cm < 0.00 

(GcV/c) [pb/(GeV/c)^] [pb/(GeV/c)^] [pb/(GcV/c)^] 



3.50 


- 4.00 


1150 ± 220 ± 210 


1470 ± 210 ± 270 


1390 ± 220 ± 260 


4.00 


- 4.50 


339 ± 14 ± 57 


00-1 1 'i c\ 1 A 

381 ± 12 ± 64 


365 ± 11 ± 62 


4.50 


- 5.00 


87.9 ± 4.3 ± 13 


105.5 ± 4.5 ± 16 


■\ A 1 A 1 T/^ 

105.4 ± 4.7 ± 16 


5.00 


- 5.50 


22.5 ± 1.9 ± 3.0 


28.2 ± 2.0 ± 3.8 


36.4 ± 2.3 ±4.9 


5.50 


- 6.50 


5.19 ± 0.48 ± 0.65 


7.41 ± 0.60 ± 0.92 


8.25 ± 0.68 ± 1.0 


6.50 


- 8.00 


0.61 ± 0.11 ±0.08 


0.97 ± 0.14 ±0.12 


1.20 ± 0.16 ±0.15 


8.00 


- 10.00 


0.020 ± 0.016 ±0.003 


0.034 ± 0.019 ±0.005 


0.038 ± 0.020 ± 0.005 






0.00 < t/cm < 0.25 


0.25 < t/cm < 0.50 


0.50 < t/cm < 0.75 


3.50 


-4.00 


1510 ± 190 ±280 


1570 ± 160 ± 290 


1120 ± 260 ± 200 


4.00 


-4.50 


370.8 ± 9.9 ± 63 


367.2 ± 9.3 ± 62 


290 ± 52 ±49 


4.50 


- 5.00 


107.2 ± 4.3 ± 16 


111.4 ± 4.2 ± 16 


87.1 ± 4.0 ± 13 


5.00 


- 5.50 


37.2 ± 2.2 ±5.0 


36.8 ± 2.1 ±4.9 


28.1 ± 2.0 ±3.8 


5.50 


- 6.50 


8.95 ± 0.65 ±1.1 


8.19 ± 0.63 ±1.0 


7.42 ± 0.61 ±0.92 


6.50 


- 8.00 


0.85 ± 0.14 ±0.11 


1.07 ± 0.15 ±0.13 


0.62 ± 0.13 ±0.08 


8.00 


- 10.00 0.055 ± 0.025 ± 0.007 


0.040 ± 0.021 ± 0.005 


0.038 ± 0.022 ± 0.005 



TABLE IV: The invariant differential cross section [Eda/dPp) per nucleon as a function of j/cm and pr for the inclusive reaction 
pBe 7X at 800 GeV/c. 





Pt 


-1.0 < < -0.75 


-0.75 < 2/cm < -0.50 


-0.50 < 2/cm < -0.25 


(GeV/c) 


[pb/(GeV/c)2] 


[pb/(GeV/c)2] 


[pb/(GeV/c)^] 


3.50 


- 4.00 


1200 ± 320 ±230 


2000 ± 320 ± 380 


2020 ± 320 ±380 


4.00 


- 4.50 


516 ± 43 ±90 


499 ± 44 ± 87 


703 ± 44 ± 120 


4.50 


- 5.00 


129.3 ± 7.4 ±20 


157.7 ± 8.2 ±24 


218.8 ± 8.9 ±33 


5.00 


- 5.50 


44.2 ± 3.8 ± 6.0 


64.6 ± 4.0 ±8.8 


78.2 ± 4.6 ± 11 


5.50 


- 6.50 


10.31 ± 0.93 ± 1.3 


18.5 ± 1.3 ±2.3 


23.3 ± 1.4 ±2.9 


6.50 


- 8.00 


1.35± 0.18±0.16 


1.93 ± 0.24 ±0.23 


3.41 ± 0.32 ±0.40 


8.00 


- 10.00 


0.170 ± 0.045 ± 0.020 


0.169 ± 0.046 ± 0.020 


0.325 ± 0.070 ± 0.038 






-0.25 < ycm < 0.00 


0.00 < Vera < 0.25 


0.25 < J/cm < 0.50 


3.50 


-4.00 


3710 ± 410 ± 700 


2680 ± 250 ± 510 


1870 ± 250 ± 350 


4.00 


-4.50 


792 ± 39 ± 140 


898 ± 38 ± 160 


763 ± 35 ± 130 


4.50 


- 5.00 


234.7 ± 8.5 ±36 


297 ± 12 ±45 


250 ± 12 ±38 


5.00 


- 5.50 


95.5 ± 4.6 ± 13 


100.5 ± 5.5 ± 14 


103.4 ± 5.6 ± 14 


5.50 


6.50 


27.2 ± 1.4 ±3.4 


30.7 ± 1.6 ±3.8 


25.1 ± 1.6 ±3.1 


6.50 


- 8.00 


3.78 ± 0.32 ±0.44 


4.19 ± 0.37 ±0.49 


3.29 ± 0.37 ±0.38 


8.00 


- 10.00 0.149 ± 0.054 ± 0.017 0.285 ± 0.076 ± 0.033 


0.412 ± 0.100 ± 0.048 
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TABLE V: The invariant differential cross section {Eda/dPp) per nucleon as a function of j/cm and px for the inclusive reaction 
TT-Be 7X at 515 GeV/c. 



Pt -0.75 < 2/cm < -0.50 -0.50 < j/cm < -0.25 -0.25 < j/cm < 0.00 
(GcV/c) [pb/(GcV/c)^] [pb/(GcV/c)^] [pb/(GcV/c)^] 



3.50 


4.00 


1250 ± 220 ±210 


1150 ± 100 ± 190 


1312 ± 88 ± 220 


4.00 


4.50 


370 ± 15 ± 58 


342 ± 11 ±53 


401 ± 11 ±63 


4.50 


- 5.00 


80.4 ± 5.0 ± 11 


124.3 ± 4.7 ± 17 


137.4 ± 5.1 ± 19 


5.00 


- 5.50 


30.7 ± 2.2 ±3.9 


48.0 ± 2.5 ±6.0 


54.6 ± 2.7 ±6.8 


5.50 


- 6.50 


8.63 ± 0.64 ± 1.0 


12.79 ± 0.79 ± 1.5 


18.51 ± 0.93 ± 2.2 


6.50 


- 8.00 


0.96 ± 0.14 ±0.11 


1 76 + n 1 9 ± n 20 


2 84 + n 24 ± n 33 


8.00 


- 10.00 


0.060 ± 0.029 ± 0.007 


0.159 ± 0.040 ±0.019 


0.294 ± 0.056 ± 0.035 






0.00 < 2/cm < 0.25 


0.25 < 2/cm < 0.50 


0.50 < 2/cm < 0.75 


3.50 


-4.00 


1670 ± 120 ± 280 


1490 ± 100 ± 250 


1240 ± 67 ± 210 


4.00 


- 4.50 


436 ± 11 ± 68 


444 ± 11 ± 69 


387± 11 ±61 


4.50 


- 5.00 


162.7 ± 5.1 ± 22 


175.2 ± 5.4 ±24 


165.5 ± 5.4 ±23 


5.00 


5.50 


64.1 ± 2.7 ±8.1 


73.2 ± 2.9 ±9.2 


69.1 ± 3.0 ±8.7 


5.50 


- 6.50 


21.73 ± 0.95 ±2.6 


23.30 ± 0.99 ±2.8 


21.9 ± 1.0 ±2.6 


6.50 


-8.00 


3.22 ± 0.25 ±0.37 


4.09 ± 0.29 ± 0.47 


3.10 ± 0.28 ±0.36 


8.00 


- 10.00 0.424 ± 0.065 ± 0.051 


0.367 ± 0.068 ±0.044 


0.237 ± 0.053 ±0.028 



TABLE VL The invariant differential cross section (^Eda/d'^p) as a function of jjcm and pr for the inclusive reaction pp — > 7X 
at 530 GeV/c. 

-0.75 < ycm < -0.50 -0.50 < y^m < -0.25 -0.25 < ycm < 0.00 



(GcV/c) [pb/(GcV/c)^] [pb/(GcV/e)^] [pb/(GcV/c)^] 



3.50 


4.00 


2210 ± 570 ±400 


1180 ± 480 ±220 


850 ± 420 ± 160 


4.00 


- 4.50 


356 ± 32 ± 60 


363 ± 28 ±61 


411 ± 25 ±69 


4.50 


- 5.00 


93± 12± 14 


112 ± 11 ± 16 


128 ± 11 ±19 


5.00 


- 5.50 


24.8 ± 5.1 ±3.3 


35.8 ± 5.0 ±4.8 


38.0 ± 5.6 ±5.1 


5.50 


- 6.50 


5.6 ± 1.2 ±0.7 


9.3 ± 1.6 ±1.2 


6.9 ± 1.6 ±0.9 


6.50 


- 8.00 


0.25 ± 0.20 ±0.03 


0.84 ± 0.32 ± 0.10 


0.87± 0.34 ± 0.11 






0.00 < ycm < 0.25 


0.25 < Vcm < 0.50 


0.50 < Vcm < 0.75 


3.50 


4.00 


1110 ± 390 ±200 


800 ± 240 ± 150 


560 ± 520 ± 100 


4.00 


- 4.50 


315 ± 23 ±53 


327 ± 22 ±55 


390 ± 160 ± 70 


4.50 


- 5.00 


110 ± 11 ±16 


103.7 ± 9.6 ±15 


72.9 ± 9.2 ±11 


5.00 


- 5.50 


36.0 ± 5.3 ±4.8 


33.0 ± 5.1 ±4.4 


24.2 ± 4.7 ±3.2 


5.50 


- 6.50 


10.4 ± 1.7 ± 1.3 


10.0 ± 1.6 ± 1.2 


4.0 ± 1.3 ±0.5 


6.50 


- 8.00 


1.01 ± 0.35 ±0.12 


1.65 ± 0.41 ±0.20 


0.60 ± 0.27 ±0.07 



TABLE VII: The invariant differential cross section {Eda/dPp) as a function of j/cm and px for the inclusive reaction pp 

at 800 GeV/c. 





Pt 


-1 < Vrrr, < "0 75 


—0.75 < j/cm < —0.50 


—0.50 < j/cm < —0.25 


(GcV/c) 


[pb/rGcV/c)2l 


[pb/(GcV/c)2l 


[pb/(GeV/c)2l 


3.50 


- 4.00 


1990 ± 870 ± 380 


1320 ± 790 ± 250 


990 ± 780 ± 190 


4.00 


- 4.50 


483 ± 72 ± 84 




744 ± 83 ± 130 


4.50 


- 5.00 


162 ± 20 ± 25 


211 ± 20 ± 32 


219 ± 21 ± 33 


5.00 


- 5.50 


44.0 ± 7.6 ± 6.0 


78.9 ± 9.6 ± 11 


84± 11 ± 11 


5.50 


- 6.50 


12.6 ± 2.3 ± 1.6 


20.8 ± 3.0 ± 2.6 


25.0 ± 3.3 ± 3.1 


6.50 


- 8.00 


1.69 ± 0.49 ± 0.20 


2.20 ± 0.64 ± 0.26 


3 02 ± 79 ± 35 


8.00 


- 10.00 


0.18 ± 0.11 ± 0.02 


45 ± 20 ± 05 


44 ± 19 ± 05 






-0.25 < t/cm < 0.00 


0.00 < 2/cm < 0.25 


0.25 < j/cm < 0.50 


3.50 


-4.00 


2290 ± 810 ±430 


2510 ± 570 ±470 


2190 ± 550 ±410 


4.00 


-4.50 


801 ± 52 ± 140 


810 ± 100 ± 140 


621 ± 79 ± 110 


4.50 


- 5.00 


208 ± 19 ±31 


330 ± 27 ± 50 


222 ± 30 ± 34 


5.00 


- 5.50 


85 ± 10 ± 12 


111± 13 ±15 


90 ± 13± 12 


5.50 


6.50 


30.4 ± 3.5 ±3.8 


31.5 ± 3.8 ±3.9 


23.7 ± 3.8 ±2.9 


6.50 


8.00 


3.74 ± 0.79 ±0.44 


5.28 ± 0.96 ±0.62 


4.18 ± 0.94 ±0.49 


8.00 


- 10.00 


0.25 ± 0.16 ±0.03 




0.77 ± 0.23 ±0.09 



TABLE VIII: The invariant differential cross section (^Eda/dPp) as a function of j/cm and px for the inclusive reaction tt p 
at 515 GeV/c. 



Pt 


-0.75 < 2/cm < -0.50 


-0.50 < 2/cm < -0.25 


-0.25 < 2/cm < 0.00 


(GeV/c) 


[pb/(GeV/c)^] 


[pb/(GeV/c)^] 


[pb/(GeV/c)2] 


4.00 


4.50 


477 ± 76 ± 74 


454 ± 79 ± 71 


566 ± 71 ± 88 


4.50 


5.00 


156 ± 26 ±21 


158 ± 27 ± 22 


189 ± 29 ± 26 


5.00 - 


5.50 


43 ± 12 ±5.0 


52 ± 14 ±6.0 


73 ± 16 ±9.0 


5.50 - 


6.50 


6.5 ± 3.5 ±0.8 


14.2 ± 4.6 ± 1.7 


20.5 ± 5.0 ±2.4 


6.50 - 


8.00 


2.1 ± 1.0 ±0.2 


2.9 ± 1.2 ±0.3 


4.0 ± 1.5 ±0.5 






0.00 < I/cm < 0.25 


0.25 < 2/cm < 0.50 


0.50 < 2/cm < 0.75 


4.00 - 


4.50 


514 ± 64 ±80 


547 ± 57 ± 85 


403 ± 66 ± 63 


4.50 


5.00 


137 ± 26± 19 


208 ± 29 ± 29 


145 ± 29 ± 20 


5.00 - 


5.50 


65 ± 14 ±8.0 


84± 16± 11 


72 ± 16 ±9.0 


5.50 - 


6.50 25.9 ± 5.3 ±3.1 


25.0 ± 5.5 ±3.0 


19.8 ± 5.5 ±2.3 


6.50 - 


8.00 


4.7 ± 1.6 ±0.5 


4.7 ± 1.6 ±0.5 


4.9 ± 1.9 ±0.6 



